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ABSTRACT 

Aims. We calculate stationary configurations of rapidly rotating compact stars in general relativity, to study the proper- 
ties of circular orbits of test particles in the equatorial plane. We search for simple, but precise, analytical formulae for 
the orbital frequency, specific angular momentum, and binding energy of a test particle that are valid for any equation 
of state and any rotation frequency of the rigidly rotating compact star, up to the mass-shedding limit. 
Methods. Numerical calculations are performed using precise 2-D codes based on multi-domain spectral methods. Models 
of rigidly rotating neutron stars and the space-time outside them are calculated for several equations of state of dense 
matter. Calculations are also performed for quark stars consisting of self-bound quark matter. 

Results. At the mass shedding limit, the rotational frequency converges to a Schwarzschildian orbital frequency at the 
equator. We show that the orbital frequency of any orbit outside equator can also be approximated by a Schwarzschildian 
formula. Using a simple approximation for the frame-dragging term, we obtain approximate expressions for the specific 
angular momentum and specific energy on the corotating circular orbits in the equatorial plane of neutron star, which 
are valid to the stellar equator. The formulae recover reference numerical values with typically 1% of accuracy for 
neutron stars with M > 0.5 Mq. These values are less precise for quark stars consisting of self-bound quark matter. 
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1. Introduction 

Circular orbits of test particles moving freely along 
geodesies in a neutron star equatorial plane repre- 
sent, to a very good approximation, orbits of gas 
elements in a low-mass X-ray binary (LMXB) thin 
accretion disk. The simplest model of these orbits 
is obtained for a Schwar zschild space-time p r oduce d 
by a static neutron star (|Shapiro fc TeukolskvL Il983h . 
However, accretion in LMXBs is associated with neu- 
tron star spin-up, and is a commonly acc epted mech- 
anism for producing millisecond pu ls ars (lAlpar et all 



Il982t iBhattacharva fe van den Heuvell . Il99lh . This sce- 
nario is corroborated by the discovery of rapid pulsa- 
tions with frequencies up to 619 Hz, in more than a 
doze n LMXBs. List of such b ursters is given in Table 
1 of iKiziltan fc Thorsettl (120091) . wh i ch alr eady needs an 
update: according to IGallowav et all (|2009D . the transient 
burster EXO 0748—676 pulsates actually at 552 Hz instead 
of 45 Hz, and a new source Swift J1749.4— 2807, pu lsating 
at 518 Hz, was reported by lAltamirano et all (|2010t ). Nine 
neutron stars in LMXBs spin at frequencies ranging from 
524 Hz to 619 Hz, and their rotation significantly affects 
the space-time and thus the orbits of test particles in their 
vicinity. 

At a given distance from the star center, a test particle 
moves at an orbital frequency / or b on a circular orbit in an 



equatorial plane, with the energy and angular momentum 
per unit rest mass E and I, respectively. Axial symmetry is 
assumed, and all quantities are defined in a chosen reference 
system, using suitable space-time coordinates. The star is 
assumed to be rigidly rotating, its gravitational mass and 
angular momentum being M and J, respectively. Analytical 
formulae expressing / or b, E, and I as functions of radial 
coordinates were obtained in the slow-rotation approxima- 
tion, keeping only linear terms in the star 's angular mo- 
mentum J, bv iKluzniak fc Wagoner! (|l985l ). Their expres- 
sions coincide, as they should, after an appropriate change 
of coordinates, with those obtained in the lowest order in 
J for the Kerr metr ic of a rotating black hole o f the same 
M and J (see, e.g.. IShapiro fe Teukolskv||1983t ). However, 
exact 2-D calculations show that the slow-rotation approx- 
imation is not valid for orbits in the vicinity of compact 
stars (neutron stars and quark stars) ro t ating at frequency 
higher than 400 Hz (I Miller et all 119981: IShibata fe Sasaki 
119981: IZdunik et all l2002t iBerti et all 120051) . 
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The present research was motivated by a specific as- 
trophysical project: modeling of the spin-up process of an 
old "radio-dead pulsar", via accretion within a LMXB , 
into a millisecond pulsar (jBhattacharva fc van den Heuvell 
I1991D . During this recycling process, the neutron-star mag- 
netic field pushes the inner radius of thin Keplerian disk to 
ro, which is usually significantly larger than the radius of 
the innermost stable circular orbit around an idealized zero 
magnetic field star. To perform efficiently a large number 
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of simulations corresponding to different astrophysical sce- 
narios, one needs reliable analytic expressions for / O rb(?"o) 
and 1{tq), valid for any rotation frequency up to the mass- 
shedding limit /k, and for ranging from the stellar ra- 
dius to a few thousand of radii. This range of tq is re- 
quired to model the complete process of the recycling of an 
old "dead" pulsar rotating at ~ 0.1 Hz, with polar mag- 
netic field ~ 10 12 G, to a millisecond pulsar rotating at 
~ 500 — 1000 Hz, and also a pulsar w ith a polar magnetic 
field - 10 8 - 10 9 G (|LorimerL 120081) . The inspiration for 
the present note came from our previous puzzling result on 
the implications of (still unconfirmed) detection of neutron 
star rotation at 1122 Hz for the equ ation of state (EOS) 
of dense matter (|Beiger et all I2007D . At a given gravita- 
tional mass M, we found remarkably that the orbital fre- 
quency in the Schwarzs child space-time produced by M, 
calculated at the actual mass-shedding equatorial circum- 
ferential radius, coincided (within a fraction of a percent) 
with a true mass-shedding frequency for this M, which is 
equal to the true orbital frequency of a test parti cle at 
the equator (|Beiger et all l2007t lHaensel et all . I2008D . The 
Schwarzschildian formula yielded true / or bj when one in- 
serted into it an actual equatorial radius. As we show ex- 
plicitly in the present note, this is due to the mutual can- 
celation of the effects of both the dragging of the inertial 
frames and the neutron star oblateness. 

In the present note, we show that using a 
Schwarzschildian formula for / rb(Po)j combined with 
simple approximations for basic quantities (metric func- 
tions, orbital velocity) one can reproduce very precisely 
(within a fraction of a percent) orbital parameters outside 
a rotating neutron star, down to the star surface, and for 
rotation frequencies up to the mass-shedding limit. We 
explain this remarkable property using existing expansions 

of /orb, E, and I i n the powers of the spec ific stellar an- 
gular momentum ([Abramowicz et all I2003D . In addition, 
using exact 2-D calculations of space-time of a rotating 
neutron star, we study how Schwarzschildian / orD at the 
equator converges, in the mass-shedding limit, to the star 
rigid-rotation frequ ency, a feature pointed out already in 
iBeieer et all (|2007t) . 

Notations, space-time coordinates, and exact expres- 
sions for the metric functions and for quantities associated 
with circular equatorial orbits are introduced in Sect. [5] 
Two well-known analytic models of equatorial orbits are 
summarized in Sect. [3] Our approximation for / or b is de- 
scribed in Sect. 21 Approximations for I and E are presented 
in Sect. [5] Approximations for / or b, E, and I are compared 
with the exact values of these quantities, and their high 
pr ecision is explained using systematic expansions derived 
bv lAbramowicz et all (|2003h . T he formula for the specifi c 
energy of test particle, given in lAbramowicz etTa l. (2003), 
contains a misprint that we correct. Finally, Sect. [B] con- 
tains a discussion of our results and conclusions. Details 
concerning equations of state, used in our calculations, are 
collected in the Appendix. 



2. Notations, metric, and orbits 

The metric for a rapidly rotating object, assuming that the 
space-time is axisymmetric, stationary, asymptotically flat, 
and free of meridional currents can be written by means of 
the 3+1 formalism of GR, in the so-called maximal-slicing 
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Fig. 1. Frequency of stellar rigid rotation (solid line) 
compared with the frequency of a test particle orbiting 
a point mass, equal to the mass of the rotating star 
(Schwarzschildian frequency), at a distance of stellar cir- 
cumferential radius (dotted line) . Both frequencies approx- 
imately coincide at the mass-shedding configuration radius 
( marked with a circle) . Stella r model based on the DH EOS 
of iDouchin fc Haensell (|200lh described in the Appendix. 

quasi-isotropic coordinates 

ds 2 = -N 2 dt 2 + A 2 {dr 2 +r 2 d9 2 ) 

+ B 2 r 2 sin 2 d{d(f> - N^dt) 2 , (1) 

where A and B are factors of the 3-metric (A = B in 
spherical symmetry), N is called the lapse, and is 
the only non-zero component of the shift vector. The in- 
tegrals of motion can be obtained by constructing the 
Lagrangian and using the variational principle (see e.g., 
Shapi ro fc Teukolskvlfl98l . There are two integrals of mo- 
tion of a test particle of rest mass m, corresponding to the t 
and 4> coordinates: E, the energy of a test particle (in units 
of mc 2 ), and I, the specific angular momentum (in units of 
mc). One can also define V, the effective potential 



V 



N 2 1 + 



I 2 



B 2 r 



+ 2N^El - l 2 



Conditions for circular orbits are 



E 2 = V 2 



and V tr = 0. 



(2) 



(3) 



The radius r C i rc of the circular orbit in the equatorial plane 
of the star (circumferential radius) is related to the coordi- 
nate radius by the relation 

Tcvrc = Br. (4) 

The angular velocity of the matter, as measured by a dis- 
tant observer, can be described by the components of the 4- 
velocity vectors u {q a BU a u^ = —I ) and equals £1 = u^/u* 
ria nt. Following iBardeenl (|l97fl . the proper velocity of 
matter in the equatorial plane is 



v = (n-N'f')- 



N 



(5) 



The specific (per unit mass) energy E of a particle and its 

specific angula r momentum / can be expressed in terms of 
the velocity v (|BardeenLll972l ): 



E 



N + A^Wc, 



(6) 
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From both expresions in Eq. , one obtains the corotating 
test particle velocity v 



v = ^Ar{lnN). r + 4r 2 (ln 7V),,.(ln B) . r + {r 2 Jj N f) 2 
+ r 2 ^^)/(2 + 2r(lnB), r ), (7) 



and, from Eq. ([5]), the orbital frequency in the equatorial 
plane 

3. Approximate analytic solutions 

3.1. Schwarzschild solution 

Simplest approximation consists of neglecting neutron star 
rotation, and dealing with the exterior space-time of sta- 
tionary, spherically-symmetric object. The metric in the 
Schwarzschild coordinates is then 



ds 2 



r 



dr 2 



l — r s /r 

+ r 2 (d0 2 + sin 2 0d(/) 2 ), (9) 



where r s — 2GM/c 2 is the Schwarzschild radius and M is 
the total gravitational mass of the object. The coordinate 
radius r is here equal to the cicumferencial radius r C i rc de- 
fined in Eq. (U). Specific energy E Schw - and specific angular 
momentum l Schw - are given by the formulae 



E 



Schw. 



(r 2 — 1.5r s r) 



1/2 ' ' 



Schw. 



I GMr 2 
r — 1.5r, 



The corresponding orbital frequency f^ b hw ' is 

J orb n _ \/ -.3 



eSchw. 



1 GM 



2tt 



(10) 



(11) 



3.2. Slow-rotation approximation 

Another popular appr oximation c an be easily obtained 
from the Kerr solution (jKerrl . I1963T) in Boyer-Lindquist co- 
ordinates; this slow-rotation approximation retains only the 
first-order terms in the angular momentum J (or equiva- 
lently, in a = J /Ale). The formulae for E sr , I sr and f£ b 
are 



rsr _ r 
J orb .1 ( 




(12) 

(13) 
(14) 



4. Approximation of orbital frequency by the 
Schwarzschildian formula 

The seemingly surprising property of angular frequency f or b 
being approximately equal to f^ b hw ' can be explained by 



studying the higher-order expansion terms of axisymmet- 
ric metrics, such as Ha rtle-Thorne metric, presented by 
HE ramowicz et all (|2003T l (hereafter referred to as AAKT). 
AAKT give the expansion to the second order in angular 
momentum and include the mass-quadrupole moment, Q, 
term, which plays a decisive role in a co rrect description of 
the orbital parameters of a test pa r ticle (Shibata fc Sasaki 
119981: lLaarakkers fc Poissonl 119991: iBerti et al l . 120051) . The 
orbital angular velocity of a corotating test particle can be 
approximated, according to AAKT, as 



CI = 



M 1 ' 2 

r 3/2 



.M3/2 
' r 3/2 



+ j 2 F 1 {r) + qF 2 (r) 



(15) 



where G = c = 1, j = a/M , q = Q/M 3 and the functions 
F\ (r) and F% (r) are defined therein. 

As an illustration, we presenlQ in Fig. 2 a compar- 
ison of the frequency fo rb KT = fl/2n, for the mass- 
shedding configura t ion c alculated for realistic EOS of 
iDouchin fc Haensel (|200ll ) (details in the Appendix), with 
Ub, fit™-, and fZ b (Eqs.lHU and M respectively). 
As expected, f^ KT reproduces the true, numerically- 
obtained values of f or b far more precisely than f^ b , es- 
pecially near the stellar surface. One has to keep in mind 
however that in Eq. (1151) we have used the mass-quadrupole 
moment obtained directly from our numerical calculations 
instead of an approximation to q, which may affect the be- 
havior of f^ r f KT . The value of f^,.f KT is also quite close to 
forb ■ We sought for an explanation of this phenomenon 
by analyzing the terms in Eq. (j!5[) : in addition to them 
being small in comparison to the leading term, the j-term 
is approximately equal, but with the opposite sign, to the 
mass-quadrupole moment g-term, so that they effectively 
cancel each other (see Fig. [5] for an typical comparison; 
the j 2 -term is much smaller than both of the other terms, 
two orders of magnitude in this example, and therefore in- 
significant). This feature is qualitatively and quantitatively 
present for different polytropic and realistic EOSs as well 
as bare strange-quark matter stars (details concerning the 
EOSs that we used can be found in the Appendix). 

The quadrupole term in Eq. (IT5|) is related to the ro- 
tational oblateness of the star and, in contrast to j-terms, 
is of Newtonian nature. One can study the maximum de- 
viation from the Schwarzschildian (and Newtonian) test 
particle orbital frequency resulting from this effect, assum- 
ing a dense matter disk (which corresponds to an "ex- 
treme oblateness") as a source of gravitational field. The 
gravitational pull is described by hypergeometric functions 
(jZdunik fc Gourgoulhonl l200ll ) that indicate the maximum 
frequency deviation A/// ~ 18% for orbital radii larger 
than that of the innermost stable circular orbit - the exis- 
tence of unstable circular orbits is a result of the oblateness 
of the gravitational field source, treated in the Newtonian 
theory. 



1 As in Bert i et "all l|2005l ) (Sect. 3.1 therein), we use exact 
numerical values of the gravitational mass M, total angular mo- 
mentum J, equatorial radius R, but also the mass-quadrupole 
moment Q of a given configuration as input to AAKT expan- 
sions to reproduce as well as possible the test particle orbital 
frequencies, angular momenta, and energies and compare them 
with their analogues from different space-times. 
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Fig. 2. Upper panel: percentage differences between ex- 
act numerical values of test particle orbital frequencies, 
/ or b, and approximations presented in the text, plotted 
against the radius of a circular orbit; A/// = (/orb — /)//, 
where / = f^ r i KT ,forb^ or fofb™'- Lower panel: values of 
various terms in Eq. (|15D versus orbital radius. For com- 
parison, relative frequency difference in the case of the 
Schwarzschildian formula is also shown (mass-shedding con- 
figuration of DH EOS, details in Appendix). 

5. Approximation of the specific orbital energy and 
angular momentum 

We approximate the specific energy and angular mo- 
mentum of a t e st pa rticle, in a manner similar to 
iShibata &; Sasaki (|1998l) . As one can see in an example 
shown in Fig. [21 the difference between the exact orbital 
frequency and the Schwarzschildian formula is smaller than 
1% even for mass-shedding, compact stellar configurations. 
Near the surface of the star, f^ b hw ' is far more accu- 
rate than the slow-rotation approximation, and as good as 

f AAKT 
J orb 

Motivated by this result, we aim to provide formulae 
that are simpler than those of AAKT and which could 
serve the ambition of being useful in practical calculations. 
To this aim, we propose to substitute the exact values 
for the test particle velocity in the Eq. ([5]) by their ap- 
proximations of the angular frequency il by 2ir fg r c b hw ' from 
Eq. (|TT|) . the azimuthal shift component by the first- 
order term in the slow-rotation approximation of a frame- 
dragging term 2GJ/r 3 c 2 , and the metric function N by its 



Schwarzschildian c^/l — r s /r: 



Vappr. — \^^fi 



Schw. 
orb 



2GJ/r 3 c 2 ) 



c^Jl - r s /r 



(16) 



where r 



for brevity. We then apply v appr , , as well as 



the approximations described above, to Eqs. ©, obtaining 
l a PP r - and E appr - . In Figs. [21 and S] we compart our approx- 
imation with the Schwarzschild, slow-rotation, and AAKT 



2 Note the difference in definitions o f the specific orbita l angu- 
lar momentum between the article of lAbramowicz et alj (2003) 
and here: l AAKT — —u^/ut = —l/E. We also notice that there is 
a misprint (?) in the Eq. (30) of AAKT. Their expression leads 
to divergence, with increasing j, of approximate E(r) from the 
exact values of this quantity. This could be repaired by replacing 
-20M 2 r in their Eq. (30) by +20M 2 r. 



(Eqs. (21) and (26) therein) formulae for I and E. For dif- 
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Fig. 3. Comparison of specific orbital angular momenta as 
a function of the radius of the circular orbit. The exact 
value is denoted by I, while l a is an approximation of it. 
The relative deviations (in percent) are Al/l = (I — 
where l a = J^ppr- Jaakt J sr QT Jschw._ Results for DH E0 S 
mass-shedding configuration and for the currently-highest 
pulsar frequency, 716 Hz (details in Appendix). 
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Fig. 4. Comparison of specific orbital energies as a func- 
tion of the radius of the circular orbit. Configurations and 
notation analogous to Fig. [31 

ferent EOSs, masses, and rates of rotation, the behaviour of 
these deviations is quantitatively and qualitatively similar. 



6. Discussion and conclusions 

We have performed calculations of stationary configura- 
tions of rotating compact stars, for a set of representa- 
tive EOSs, both polytropic and realistic ones, as well as 
for the bag models describing the EOSs of hypo thetical 
bare strange quark stars (as in Hacnscl et al.ll2009r ). In all 
these cases, we observed quantitatively and qualitatively 
similar behaviour of the orbital frequency of a test particle 
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moving on a circular orbit in the equatorial plane. In all 
cases, the value of the orbital frequency was close to that 
of a test particle in a Schwarzschildian space-time around 
a point mass (equal to that of the actual compact star), 
at the properly defined radius. This result is valid for any 
EOS and for any rotation frequency of the compact star, 
up to the mass-shedding limit. 
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Fig. 5. Comparison of orbital frequencies for stellar con- 
figurations of DH EOS and bare strange quark star EOS, 
rotating at 716 Hz. Notations as in Fig. [21 details in 
Appendix. 



derivatives of metric functions (see, e.g.. lCook et alJll994fl . 




Fig. 6. Comparison of I, E, and f orb for 2.2 M Q SQM1 EOS 
star rotating at 1300 Hz. Notation as in previous figures. 



We note that the approximations of I and E, proposed 
in Sect. [5J reproduce true (numerical) values within about 
one per cent (being the least accurate in the vicinity of 
the stellar surface). Their accuracy near the stellar surface 
decreases to a few percent for rapidly rotating stars with 
extreme compactnesses (2GM / Rc 2 > 0.5), especially in the 
case of rapidly rotating and oblate bare strange-quark mat- 
ter stars with large quadrupole moments, but as is shown in 
Fig. [5] the approximations are acceptable for the currently 
highest observed spin rate of 716 Hz. For comparison, Figs. 
[3]and[2]show the difference in the precision for the DH EOS 
in the case of the mass-shedding limit and the rotation fre- 
quency 716 Hz. 

Overall, our approximations for I and E are much bet- 
ter than the slow-rotation approximation, and one does 
not need to compute the mass-quadrupole moment of 
the star, the knowledge of which is otherwise c rucial 
in the existing sy s temat ic expansions (Be rti et al.l . [2005; 
lAbramowicz et all 120031) . In cases of larger compactness 
and oblateness, as well as sub-millisecond rotation periods, 
the approximation from Sect. [5] is in closer agreement with 
the results of numerical simulations than the formulae of 
lAbramowicz et al.l (|2003l ): in Fig. El we present an extreme 
example of a bare quark star spinning at 1300Hz. 

The analytic formulae for / rb( r )i K r )i an d E(r) give a 
very good approximation (typically within a one per cent) 
of exact values for neutron stars with astrophysically rele- 
vant masses, M > 0.5 M Q . They are valid for rigidly rotat- 
ing neutron stars, for rotation rates up to the mass shedding 
limit, and down to the equator. However, they cannot be 
used to calculate the properties of the innermost stable cir- 
cular orbit, because this calculation involves second radial 
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Appendix: Numerical implementation and the EOSs 
chosen for calculations 

The calculations have been performed using the rotst ar code f rom t he 
LORENE librarjB We use EOSs employed in lHaensel et al.l ||2009I1 : 

1. Realistic microphysical EOSs of dense matter. Unless marked 
otherwise, the figures e mploy the results for DH EOS of 
iDouchin fc Hacnscl (2001), for a stellar configuration correspond- 
ing to a non-rotating star of mass M = 1.43 Mq (central baryon 
density n,;, c = 0.55 fm -3 ). The mass-shedding configuration has 
a gravitational mass of 1.8 Mq, whereas that rotating at 716 Hz 
has a mass of 1.52 Mq . 

2. Bag models of bare strange quark stars. Figs. [5] and [6] show re- 
sults for SQM1 EOS. The non-rotating counterpart to the stellar 
configuration presented in Fig. [5] has the gravitational mass of 
1.43 Mq. The Fig. [6] configuration has a central baryon den- 
sity of rij, c = 0.63 fm — 3 , non-rotating gravitational mass of 
M = 1.7 Mq, and compactness 2GM/Rc 2 = 0.47. 

3. Relativistic po lytropes. We considered the range 7 = 1.5 — 3 
llTooped.ll965n . 
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